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METHOD AND APPARATUS OF NOISE VARIANCE ESTIMATION FOR 

USE IN WIRELESS COMMUNICATION SYSTEMS 



Field of the Invention 

The present invention relates generally to a method and apparatus of 
noise variance estimation for use in wireless communication systems, a nd 
more particularly, to a method and apparatus of noise variance estimation by 
exploiting the training sequence. 

Background of the Invention 

CDMA (Code Division Multiple Access) is a new wireless 
communication technology developed after FDMA (Frequency Di vision 
Multiple Access) and TDMA (Time Division Multiple Access). In CDMA 
wireless communication, different UEs (user equipments) are allocated with 
different orthogonal spreading codes, and signals spread by different UEs 
with different spreading codes ca n be transferred on the same frequency 
band. 

A CDMA downlink transmission model is put forward in the paper 
entitled "Data Detection Algorithms Specially Designed For The Downlink of 
CDMA Mobile Radio Systems", VTC, 1997, by A. Klein, as shown in Fig.1 . In 
order to transmit signal vectors d (1) , ...,d (k) , d cK) (wherein d (k) (k=1 ...K) 
is composed of N complex components) to UE1, UEk, ...UEK 

respectively, base station 200 first spreads signal vectors d (1) , . . . , d (k) , . . . , d 
(K) by exploiting spreading codes c d (1) , c d (k) , c d <K) allocated to 
UE1, UEk, UEK, then combines the spread signal vectors into signal 
vector s d and transmits it to each corresponding UE 220 via the same 
channel 210. Assumed that signal vector s d reaches UEac (k=1 ...K) through 
multiple propagation paths and the CIR (channel impulse response) of each 
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propagation channel is h d a> (i=1 » 2> ... ) , signal vector e d received by 
UEk can be expressed by equation (1) as follows: 

e d (K) =H/ K) C d d+n d (K) = H d (K) s d +/7 d (K) (1) 

wherein H 6 iKy is the CI R matrix constructed with the CI R h 6 co M (i=1- 
2> ...) of each propagation channel, C d is the spreading code matrix 
constructed with spreading codes c d (1> , c d (k) , c d (K) (as to the 
construction methods of H d (/c) and C d , referring to the above paper by A. 
Klein), d=(c[ (1)T , d (k)T , c( <K>T ) T , [.] T represents matrix transposition, 
s d represents the obtained signal vector after d is spread and combined, s d 
= C d d, and n d (K) is the noise vector. 

Equation (1) indicates that the received signal vector e d °° contains 
UEk's desired signal vector d (K) , as well as signal vectors sent to other UEs 
by the base station and the noise vector. 

To help UE/c to obtain its desired signal vector d CK) from the received 

signal vector e d (/ ° with the minimum error, many method for signal reception 
have been presented, which can be referred to "Iterative Multiuser Receiver 
/Decoders With Enhanced variance Estimation", VTC, 1999, by Kimmo 
Kettunen, and "Zero Forcing an Mininum Mean -Square-Error Equalization for 
Multiuser Detection in Code -Division multiple -access channels", IEEE 
Transactions on Vehicular Technology, vol.45, pp.276 -287, May 1996, by A. 
Klein. But these methods for signal reception all rely heavily on the channel 
information, or namely noise variance, to obtain the desired signal vector 
from the received signal vector, and thus the noise variance needs to be 
computed precisely to obtain the desired signal vector with minimum error. 

To get an accurate noise variance, various noise estimation meth ods 
have been put forward. For example, a conventional variance estimation 
technique for use in AWGN channel is raised in "A novel variance estimator 
for turbo-code decoding", Proc. Of ITC , 97, pp173-178, April 1997, by M. 
Reed and J. Asenstorfer; a Rake technique for alleviating multipath 
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interference is put forward in US.PAT US2002201 10199, entitled "Method for 
Noise Energy Estimation in TDMA Systems Additionally, there are some 
noise estimation methods in which noise variance is computed by convolving 
the training sequence. These noise estimation methods can meet the 
precision requirement of 2G wireless communication systems. 

But in 3G wireless communication systems, more accurate noise 
variance is needed for signal reception, for example, the key technologies of 
multiuser detection and turbo -code both have high requirement for accurate 
noise variance. Current noise estimation methods can 't satisfy the precision 
requirement for noise variance of 3G wireless communication systems. 

Summary of the Invention 

An object of the present invention is to provide a method and apparatus 
of noise variance estimation for use in wireless communication systems, in 
which the training sequence is exploited to compute noise variance to obtain 
more accurate noise variance. 

A method of noise variance estimation is proposed in the present 
invention for use in wireless communication systems, comprising steps of: 
receiving a signal vector containing training sequence and noise vector 
transmitted via at least one propagation path f rom the base station; 
estimating, according to the signal vector, the channel impulse response of 
each propagation path to construct a channel impulse response matrix; 
calculating the noise variance of the signal vector according to the channel 
impulse response matrix and the signal vector if the channel impulse 
response remains primarily unchanged during the special time duration of 
the training sequence. 

Brief Description of the Drawings 

Fig.1 illustrates conventional CDMA downlink transmission model; 
Fig. 2 is a flow chart illustrating the noise variance estimation method in 
the present invention; 
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Fig. 3 is a block diagram illustrating the UE equipped with the noise 
variance estimation apparatus in an embodiment of the present invention; 

Fig.4 is a block diagram illustrating the noise variance estimation 
apparatus in an embodiment of the present invention. 

5 Detailed Description of the Invention 

TD-SCDMA will be exemplified in the following to describe an 
embodiment of the present invention in detail. 

In TD-SCDMA, the base station transmits signal vector to each UE in 
corresponding timeslot. According to the timeslot format of TD -SCDMA, the 

10 signal vector sent to each UE by the base station in corresponding timeslot is 

composed of the training sequence and the s pread user signal. 

With regard to the UEs allocated in the same timeslot, the base station 
first combines the signal vectors to be transmitted to each UE into a 
combined signal vector, and then transmits this combined signal vector in the 

15 timeslot to each UE. Said combined signal vector is also composed of user 

signal and training sequence, wherein the user signal in the combined signal 
vector is obtained by combining the spread user signal in the signal vector to 
be transmitted to each UE, and the training sequence in the combined signal 
vector is obtained by combing the training sequence in the signal vector to 

20 be transmitted to each UE. 

The training sequence allocated to each UE in a cell is obtained through 
performing different shift operation on the sam e basic training sequence, so 
the training sequence of the combined signal vector can be considered as 
the basic training sequence. Each UE has acquired the basic training 

25 sequence used by its cell during cell search procedure, so the training 

sequence sent by the base station in the timeslot is known beforehand to 
each UE. 

Let's suppose that the training sequence included in the signal vector 
sent by the base station in a timeslot reaches a UE through at least one 
30 propagation path, the signal vector receiv ed by the UE in the timeslot is r, 
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composed of said training sequence and noise vector n, and the known 
value of said training sequence is s. According to equation (1), signal vector r 
can be expressed as follows: 

r=Hs + n (2) 

wherein H is the CIR matrix constructed by the CIR of each propagation 
path between the UE and the base station. 

According to the channel estimation method as described in "Low Cost 
Channel Estimation in the uplin k receiver of CDMA mobile radio systems", 
Frequenz, vol.47, pp.292 -298, Nov./Dec. 1993, by B. Steiner and P.W.Baier, 

the maximum likelihood estimated value s of the training sequence 
included in signal vector r can be expressed as fol lows: 

s = (BP H) 1 U H r = s + (BfH) 1 H"n =s+n ( 3 ) 

wherein superscript H represents complex conjugate transposition. 

From equation (3), according to the known value s of the training 
sequence contained in signal vector r, the estimated value n ' of noise vector 
n can be given by: 

n =s-s=(lf H) _1 H"n (4) 

With the covariance matrix being: 

E{ri ri H } = E{(H H B.y l U H n ■ H(H^H) 1 } 

= (H^HT 1 H*£(nn" )H(H^H)" 1 } 

=c^(H // Hf 1 (5) 

wherein E{.} denotes expectation operation. By carrying out the 
operation of matrix trace between the two sides of above equation (5), it is 
easy to come down to following formulation computing the average variance 
o 2 of the estimated value n of the noise vector n: 

n 

oFj =a 2 'trace{(U H Jiy 1 }/N (6) 

wherein N is the chip duration of the training sequence , operator trace 
( ) means the computation of a matrix trace, <y 2 is the noise variance of 
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the signal vector r. 

If a w 2 is computed with conventional met hods, it will be very 

complicated. In fact, the computation of variance can be approximated 

by calculating the mean squared value of all elements about the estimated 
5 value n of the noise vector n located in one training sequence time duration 

if the channel could be regarded as constant at that time. The noise variance 
a 2 of the signal vector r can now be deduced as : 

a 2 « (n '* n ' ) / trace {(H^H)" 1 } (7) 

To further improve the estimation performance, we can sum and 
10 average the noise variance o 2 calculated from equation (7) in the timeslot 

and the noise variance o 2 calculated from equation (7) in each previous 
timeslot, and take the mean of different af as the noise variance a 2 of 

signal vector r in the timeslot. 

The above section describes the principle of computing noise variance 
15 by exploiting training sequence in the present invention. 

The following section will describe the proposed noise variance 
estimation method in detail, in conjunction with Fig.2. 

First, the UE receives a signal vector containing training sequence and 
20 noise vector in a timeslot transferred through at least one propagation path 

from the base station (step S10). 

Secondly, the UE estimates the CIR of each propagation path according 
to the received signal vector, and constructs a CIR matrix H by using the 
estimated CIR of each propagation path (step S20). 

25 Thirdly, the UE estimates the maximum likelihood es timated value s of 

the training sequence included in said signal vector using equation (3), 
according to said signal vector and said CIR matrix (step S30). 

Fourthly, the UE computes the estimated value n ' of the noise vector 
contained in said signal vector by using equation (4), according to the MLE 
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(maximum likelihood estimate) value s of the training sequence included in 
said signal vector and the known value of the training sequence (step S40). 
Wherein, the known value of the training sequence contained in said signal 
vector is acquired by the UE in cell search procedure. 

Fifthly, the UE computes the noise variance a 2 of said signal vector by 
using equation (7), according to the estimated value n ' of the noise vector 
contained in said signal vector and said CIR matrix H (step S50). Wherein 
first the power p n 2 of n' can be computed according to equation p n 2 =(n') H (n'); 
then the trace cf of matrix ( ( H H H ) can be computed, that is cf= trace (( H H H ) 

" 1 ) ; lastly, the noise variance o 2 can be computed according to equation a 2 = 
p n 2 /cf, that is , equation (7). 

Lastly, the UE sums and averages the noise variance a 2 calculated 
from equation (7) in the timeslot and the noise variance o 2 calculated from 
equation (7) in each previous timeslot, and takes the mean of different af 
as the noise variance a 2 of signal vector r in the timeslot (step S60). 

A detailed description will be given below to the proposed noise 
variance estimation apparatus, in conjunction with Fig. 3 and Fig.4. 

Fig. 3 is a block diagram illustrating the UE equipped with the proposed 
noise variance estimation apparatus. As Fig. 3 shows, in cell search 
procedure before the UE communicates with the base station, cell searching 
means 40 acquires the basic training sequence s used by the cell where the 
UE is camping. When the UE communicates with the base station, the 
antenna of the UE first sends the sign al vector Rx received in a timeslot to 
multiplier 10, and multiplier 10 multiplies the received signal vector Rx by the 
RF carrier generated by VCO 20, to convert signal vector Rx into baseband 
signal vector. Then, ADC 30 converts the baseband signal vecto r outputted 
from multiplier 10 into digital baseband signal vector r. Afterwards, cell 
searching means 40 synchronizes the digital baseband signal vector r 
outputted from ADC 30, and channel estimating means 50 computes the CIR 
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of each propagation channel for the synchronized digital baseband signal 
vector r by using conventional channel estimation methods, and constructs 
CIR matrix with the computed CIR of each propagation path. Next, noise 
variance estimating means 60 computes the noise variance of the di gital 
baseband signal vector r according to the CIR matrix computed by channel 
estimating means 50, the digital baseband signal vector r outputted by ADC 
30 and the basic training sequence s acquired by cell searching means 40. 
Finally, data detecting mean s 70 acquires the desired user signal from the 
digital baseband signal vector r according to the noise variance computed by 
noise variance estimating means 60, by using conventional data detection 
methods, such as multiuser detection method, turbo -code decoding and etc. 

Fig.4 is a block diagram illustrating noise variance estimating means 60. 
Referring to Fig.4, noise variance estimating means 60 comprises: 

equalizing means 601, for estimating the MLE value s of the training 
sequence contained in said digital baseband signal vector r according to the 
CIR matrix H computed by channel estimating means 50 and the digital 
baseband signal vector r outputted by ADC 30, by using equation (3); 

noise estimating means 602, for calculating the est imated value n' of 
the noise vector contained in said digital baseband signal vector r according 

to the MLE value s of the training sequence contained in said digital 
baseband signal vector r computed by equalizing means 601, and the basic 
training sequence s (or namely the known value of the training sequence 
contained in said digital baseband signal vector r), by using equation (4); 

noise power calculating means 603, for calculating the power p n 2 of the 
estimated value n ' of said noise vector acco rding to the estimated value n ' of 
the noise vector contained in said digital baseband signal vector r computed 
by noise estimating means 602, by using equation p n 2 =(n')" (n'); 

equalization revising means 604, for computing the trace cf of matrix 
((H"H) - 1 ) , that is cf= trace ((H H H) " 1 ) ; 

noise power revising means 605, for calculating the noise variance o 2 
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according to the power p n 2 of the estimated value n ' of said noise vector 
calculated by noise power computing means 603 and the trace cf computed 
by equalization revising means 604, by using equation a 2 = p n 2 /cf. 

Beneficial Results of the Invention 

5 As described above, in the proposed noise variance estimation method 

and apparatus for use in wireless communication systems, training 
sequence is used to compute the noise variance, so the computed noise 
variance can meet the requirement for higher accuracy. 

It is to be understood by those skilled in the art that the method and 
10 apparatus of noise variance estimation for use in wireless communication 

systems as disclosed in this invention can be modified considerably without 
departing from the spirit and scope of the invention as defined by the 
appended claims . 



